Solutions of Fisher and Burgers' equations with finite transport memory 
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The Fisher and Burgers' equations with finite memory transport, describing reaction-diffusion and 
convection-diffusion processes, respectively have recently attracted a lot of attention in the context 
of chemical kinetics, mathematical biology and turbulence. We show here that they admit exact 
solutions. While the speed of the traveling wavefront is dependent on the relaxation time in Fisher 
equation, memory effects significantly smoothen out the shock wave nature of Burgers' solution, 
without making any influence on the corresponding wave speed. We numerically analyze the ansatz 
for the exact solution and show that for the reaction-diffusion system the strength of the reaction 
term must be moderate enough not to exceed a critical limit to allow travelling wave solution to 
exist for appreciable finite memory effect. 

PACS numbers: 87.10.+e, 87.15.Vv, 87.23.Cc, 05.45.-a 



Introduction: A number of nonlinear phenomena in phys- 
ical (!]] , chemical [|| and biological processes || ^ are de- 
scribed by the interplay of reaction and diffusion or by the 
interaction between convection and diffusion. The well- 
known partial differential equations which govern a wide 
variety of them are Fisher || and Burgers' Q equations, 
respectively. While the Fisher equation describes the dy- 
namics of a field variable subject to spatial diffusion and 
logistic growth, Burgers' equation provides the simplest 
nonlinear model for turbulence. Since spatial diffusion 
is common to all these processes, Fick's law forms the 
key element in the description of transport. This de- 
scription however, gets significantly modified when the 
memory effects are taken into account, i.e., when the 
dispersal of the particles are not mutually independent. 
This implies that the correlation between the successive 
movement of the diffusing particles may be understood 
as a delay in the flux for a given concentration gradient. 
Over the last several years the analysis of memory effects 
in diffusive processes have attracted a lot of attention 
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in chemical kinetics, mathematical biology and allied ar- 
eas. The focal theme lies in the interesting traveling 
wave front solutions which have been studied extensively 
by several authors under various approximations. The 
object of the present paper is to show that the Fisher 
equation and the Burgers' equation with finite memory 
transport admit exact solutions. We numerically clarify 
the nature of the ansatz wherever necessary and analyze 
the physical implications of the solutions modified by re- 
laxation effects and the related issues. 
Fisher and Burgers ' equations with finite memory trans- 
port: The starting point of our analysis is the Cattaneo's 
modification IM of Fick's law in the form: 



which takes care of adjustment of a concentration gradi- 
ent at time t with a flux J(x, i + r) at a later time (i + r), 
t being the delay time of the particles in adopting one 
definite direction of propagation. Here u(x,t) denotes 
the field variable and D is the diffusion coefficient of the 
particles. 

The population balance equation for the particles on 
the other hand takes into account of the conservation 
equation supplemented by a source function kf(u) for 
the particles in the form 



du{x, t) 
dt 



0J_ 

dx 



kf(u) 



(2) 



The Fisher source function f(u) = u(l — u/K) has been 
the subject wide interest in various context. Here the 
first term in f(u) signifies the linear growth followed by 
a nonlinear decay due to the second one; k and K being 
the growth rate constant of the population and carrying 
capacity of the environment, respectively, fn what fol- 
lows we shall be considered with two specific cases of the 
flux-gradient relation(l) for Fisher and Burgers' problem. 
A. The Fisher equation with nonlinear damping and finite 
transport memory: We start with an expansion of J in 
Eq.(l) |p5| upto first order in r to obtain 



(3) 



Here u(x,t) represents the density function. Differenti- 
ating (3) with respect to x and differentiating (2) with 
respect to t and eliminating J from the resulting equa- 
tions one has 



J(x,t + r) = -D 



du(x, i) 
dx 



in 



£ + p -vol* 



where we have used the following abbreviations 
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/3=l/r and w 2 = (3D (5) 
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Eq. (4) , a hyperbolic reaction - diffusion equation is a gen- 
eralization of Fisher equation for finite memory transport 
and nonlinear damping, ft reduces to standard Fisher 
equation for t=0. Over the years the equation has drawn 
wide interest in the context of trave ling wave solutions in 
various problems g 1 1 0, @ MM III HM , 
g §0|, g|, g [2§. For example, Gallay and Raugel || |9|] 
have studied the propagation of front solution without 
the nonlinear term kf'{u). Horsthemke has discussed 
some related issues in the problem of transport-driven 
instabilities |l6fl . 

We now look for the traveling wave solutions of Eq.(4) 
of the form u (x, t) — KU(x—ct) = KU(z) with z = x—ct, 
where c > is the speed of the nonlinear wave (which, in 
general, is different from the linear wave w dictated by 
the medium subject to the boundary condition: 



[/(-co) = 1 and U(+oo) = 



(6) 



Eq. (4) therefore after some algebra assumes the following 
form 



d 2 U 
dz 2 



a^U „ A TT dU 

+ c(n - A)]— 2AcU — 

oz oz 



+nAmU(l -U)=0 = L(U) (say) 



where 



m = w 2 — c 2 and n — (3/m .and A — k/m 
Following Murray Q we now introduce the ansatz, 

U(z) = ' 



[1 + a exp(fez)]' 



(7) 



(8) 



(9) 



as a solution to Eq.(7), where a, 6 and s are positive 
constants to be determined. Using (9) in (7) we obtain 
after some algebra 

[s(s + l)a 2 b 2 + nAma 2 - s[c(n - A)]a 2 b 
-sa 2 b 2 ]e 2bz + [2aAmn - sab 2 - s[c(n - A)]ab] e bz 
+nAm - 2Acsabe bz (l + ae bz y s+1 
~nAm{l + ae bz )- s+2 = = L{U) . (10) 

Now for L(U) — for all z, the coefficients of e°, e bz , 
e 2bz and e 3bz within the curly brackets must vanish iden- 
tically. This implies that s—0, 1 or 2. s=0 is not a possi- 
ble solution since s is a positive constant by our starting 
assumption. For s=l the coefficients of e bz and e 2bz of 
Eq.(10) yield the following relations, 



s(.s + l)b 2 + nAm - s[c(n - A))b - sb 2 



(11) 



nAm- sb 2 - s[c(n- A)]b - 2Acsb = (12) 



which can be solved to give 6=0 and b = —2Acs/ (s + 1). 

Again since by initial assumption b is a positive con- 
stant, both the values of b are unacceptable and s=l is 
not a correct choice. 

For s=2 Eq.(10) reduces to a form in which the co- 
efficient of e bz , e 2bz and e 3bz must satisfy the following 
relations 



s(s + l)b 2 + 3nAm - 2s[c(n - A)]b - 2sb 2 = , (13) 

s(s + l)b 2 + nAm - s[c(n - A)]b - sb 2 = and (14) 

2nAm - sb 2 - s[c(n - A)]b - 2Acsb = . (15) 
From Eq.(13)-(15) we obtain 



b 2 = nAm/[s(s + 1)] 



(16) 



and putting n = /3/m, A — k/m from (8) and s=2 in 
(16) we have 



b 2 = 0k /(6m) 



(17) 



Making use of (17) in (14) we obtain b in terms of c as 
follows. 



c = 5k/3/[6b(/3-k)] 



6 = 5/ 



-<B> 



(18) 



(19) 



The exact speed c of the traveling wave can be calculated 
from (17) and (19) using m = w 2 — c 2 as 



[i + i(y-i/y) 2 ] 1/2 



(20) 



with y — y £ . It may be noted that the exact value of c 
thus derived is always greater than c m i n where 



(21) 



[l + \(y~l/yy 



,1/2 



Again in the diffusive limit, i.e., 1//3 — > or 1/y — > 
the expression (20) results in exact Fisher value of c as 
c = 5VkD/V6. It is necessary to stress that this is not 
the speed selected by the front (c = 2\JkD), but it yields 
2.04Vfc-D which is very close to the selected value [[|. 
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FIG. 1: A plot of trvelling wave solutions for different values 
of relaxation time r(= ^) for k — 0.6 and D = 1.0. The solid 
lines are due to numerical simulation of Eq.(4) and the dotted 
lines are the analytic results (22). (a) r = 0.2, (b) r = 0.4 (c) 
r = 0.6 and (d) r = 0.0 (Units are arbitrary). 



Having determined b and s one can write down the ex- 
act form of the traveling wave solution (9) for the problem 



U(z) = 1/ < 1 + oexp 



(22) 



Furthermore a can be determined from the usual condi- 
tion U{z) = 1/2 for z = 0. This results in a = (-y/2 - 1). 
The exact solution of the Fisher equation can be recov- 
ered from (22) in the limit 1//3 — > (i.e., 1/y — ► ) using 



Fisher value of c = 5VkD/\/6. This is given by 



U(z) = 



l + (V2-l)exp( A /-^ 



(23) 



We thus observe that the effect of memory or finite 
relaxation time enters into the dynamics of the reaction- 
diffusion system through its influence on the speed of the 
traveling wave front c. We emphasize here that for A = 
Eq.(22) does not give the solution selected by the front 
but is much steeper although the speed is very close to 
the selected one. 

It is pertinent to point out that although exact the 
travelling wave solution (22) does not exhaust the possi- 
bility of other solutions. This was noted earlier by Mur- 
ray Q in the context of fisher equation without memory 
effect which is a parabolic differential equation. For an 
understanding of the nature of the travelling wave solu- 
tion where (3 = (1/t) is a new element of the present 
theory, we carry out a numerical investigation of Eq. (4) 
using finite difference method to solve the boundary value 
problem. The initial condition to integrate numerically 
is that the front is at rest at t = 0. We fix the value 
of diffusion coefficient D = 1.0 for the entire treatment. 
In order to allow the variation of r for a fixed value of 
k, we have kept k at 0.6. For a higher value of k, i.e, 
where the reaction term dominates r must be chosen ap- 
propriately over a range to generate numerically stable 
travelling wave front solution. The interplay of (3 and k 
will be considered in more detail in the later part of this 
section. 

In Fig.l we compare the analytical (dotted) and the 
numerical (solid) solutions corresponding to (22) and (4), 
respectively for different values of r. From our analysis 
it is apparent that they agree fairly well for r roughly in 
the range between 0.1 and 0.5. In Fig. 1(d) we present 
the result for r — 0, which corresponds to the typical 
Fisher case. The analytical curve is marginally steeper 
than numerical one. In Fig. 2 we compare the speed of 
the travelling wave front computed numerically from (4) 
with that obtained analytically following (20) for several 
values of r. It follows that they agree reasonably well 
when r < k, i.e, in the range 0.1-0.5. As r approaches 
zero the analytical value of c becomes lower than the nu- 
merical one. This implies that the analytical wave front 
although moves slower is steeper than the numercal one 
since steepness goes as ~ - as noted by Murray in his ear- 
lier analysis. For higher values of r the disagreement be- 
tween analytical and numerical values of c grows rapidly. 

The above analysis suggests that there is a strong in- 
terplay of k and r (or (3) in the dynamics so far as the 
form and stability of the travelling wave front solution is 
concerned. To explore this aspect more clearly we now 
carry out an asymptotic analysis of the problem. To 
this end we return to Eq.(7) subject to boundary condi- 
tion (6). Following Murray we choose the perturbation 
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FIG. 2: A plot of the speed c of the travelling wave front solu- 
tion vs relaxation time r (analytical, dotted line; numerical, 
solid line) for D = 1.0, k = 0.6 (Units are arbitrary). 



Eq.(29) precludes the possibility of an explicit solution 
for <7o(£)- Depending on (5 and k we therefore consider 
three different cases; 

(i) (3 » k (or t « k): We have from (30) I = and 
(29) reduces to 

<7„(0 = (1 + exp(fcO)- 1 
or, U(z) = (1 + exp(kz/c))- 1 + O(e) (31) 

This is the standard assymptotic solution for U(z) for 
which the effect of memory is negligible. 

(ii) /3 w k (i,e r » k): We obtain similarly from (29) 
and (30) 



parameter e = 1/c 2 and look for the asymptotic solu- 
tion for < e << 1 by introducing a change of variable 
£ = - = e 1 / 2 ^ and t/(z) = ,g(0- With these transforma- 
tions Eq.(7) and (6) therefore reduces to 



2 +(n 



and 



A + 2Ag)^+mnAg(l-g) = (24) 



= i 



U(z) =11- 



exp(/c£/2) 



exp((fcz)/(2c)) 



O(e) (32) 



When both [3 and k are small compared to 1 and the 
exponential term in (32) is small it is easy to put ap- 
proximately the 0(1) term in the form of (31) as 



ff(-C») = 1 



g(+co) = 



(25) 



respectively, e in the highest derivative in Eq.(24) iden- 
tifies it as a singular perturbation problem. 

Making use of a regular perturbation series in e 



^^(i+gg^M)" 1 (33) 



(iii) (i «k (i,e t » k): We obtain 



(26) 



in (24) we obtain after equating the appropriate powers 
of e 



(n-A + 2Ag )^ = -mnAg (l - g ) ;0(1). (27) 
and 



(n-A+2A9o)^+^+2A9i^+mnA 3l (l-2 5o ) = 0; 

(28) 

The lowest order equation (27) when integrated yields 



go(0 



1 ± y/l- exp(/30 



O(e) 



(34) 



In 



(do) 13 - 



(i-sor 



(29) 



where Z is a constant of integration. Since we are inter- 
ested in the solution in the vicinity of z = 0, i.e, £ = 
for which we put <?o(£) = 1/2, we obtain 



(I)/3-fc ' 

l (IF 7 * 



■In 



(30) 



The form of this solution is generically different from 
those of (32) and (31) since it is independent of k. 

The three cases discussed above clearly shows that 
monotonic solutions satisfying U{— oo) = 1 and U(oo) = 
for finite wave speed (c > c m i n ) exist for the cases (i) 
and (ii), i,e , when r is short but finite; r < k. The as- 
sertion of this asymptotic analysis is in clear agreement 
with our numerical simulation and our choice of a smaller 
value of k as discussed earlier. 
O(e) The aforesaid analysis clearly demonstrates that al- 
though the nature of the partial differential equation 
changes from parabolic to hyperbolic type due to the 
inclusion of relaxation time, the Fisher equation can be 
solved by Murray's ansatz || to derive the exact wave 
speed and the traveling wave front solution for a suitable 
range of relaxation time r allowed by the strength of the 
reaction term. A compromise between the exact and the 
numerical solution can be obtained for relatively small 
reaction terms. The method can be extended further to 
study other density dependent diffusive processes. 
Burgers' equation with finite memory transport: The 
Burgers' equation B is a simple model of turbulence 
which illustrates an interaction between convection and 
diffusion. The convection incorporates nonlinearity in 
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the dynamics. To include finite memory effect we pro- 
ceed as follows: 

We start with the following functional relation between 
flux J{x, t+r) at a time t + r and the field variable u(x, t) 
and its gradient term at an earlier time t; 



J(x,t + r) = -u 2 (x,t) - - : - 



du(x, t) 
dx 



(35) 



where 7 is a constant. Expanding J again upto first 
order in r and differentiating the resulting equation with 
respect to x followed by differentiation of Eq.(2) for k = 
(i.e., in the absence of any source term) with respect to 
time t and elimination of J as done in the last section we 
obtain 



du 
~dt 



du 
dx 



d 2 u 



7 



dx 2 



(36) 



For r 



Eq.(36) assumes the form of classical Burgers' 



equation |l], |(| when u(x, t) and 7 are identified as the 
velocity field and kinematic viscosity, respectively. 

We now seek a traveling wave solution of the Burgers' 
equation with memory (36) in the form, U(z) — u(x — 
ct), z — x — ct, where c is again the wave speed to be 
determined. This results in the following equation: 



d 2 u TT du du n 

^ + f/ ^- c ^ = ° (37) 



where [3 = 1/t. 

We now impose the bound condition on U(z) that it 
asymptotically tends to constant values u\ as z — > —00 
and U2 as z — ► +00 and u\ > u<i- 

A direct integration of (37) yields 



dU 



1 



dz 2(4+7) 



(U 2 - 2cU - 2A) 



(38) 



where A is the integration constant. If u\ and u-i are the 
roots of the quadratic equation U 2 — 2cU — 2A — 0, then 
the wave speed c and the constant A can be obtained as 



c = Ul „ — - and A 



1 



U1U2 



2 2 
Eq.(38) can then be rewritten in the form 



(39) 



where S is given by 



5 = 



Ml - U2 



(42) 



The above analysis shows that the shape of the wave 
form is not only affected by kinematic viscosity 7 by also 
by an additional contribution c 2 / (3 due to finite relax- 
ation time t(= 1//3) such that (c 2 / j3) + 7 behaves as 
effective kinematic viscosity. It is thus apparent that the 
balance between the steepening effect of the convection 
as well as smoothening effect due to kinematic viscosity 
is enhanced by the presence of the wave speed dependent 
term c 2 //3. Thus although wave speed c [(ui + «2)/2] it- 
self remain unaffected by the finite memory effect in con- 
trast to our earlier case of Fisher equation, transmission 
layer thickness '8' - which is a measure of shock thick- 
ness, increases for higher speed c and relaxation time 
t. This implies that as the wave moves faster the shock 
smoothens out more and more so that the speed depen- 
dence of thickness S makes the dynamics self-regulating 
in the problem of interaction between convection and dif- 
fusion. 

Conclusions: The existence of relaxation or delay 
time is an important feature in reaction-diffusion and 
convection-diffusion systems. In this paper we have 
shown that two prototypical representatives of these sys- 
tems a generalized Fisher equation and Burgers' equation 
can be solved exactly for finite arbitrary delay time using 
conventional methods. While the wave speed is signifi- 
cantly modified in the Fisher problem for finite memory 
transport, speed of the traveling wave in the correspond- 
ing Burgers' problem remains unaffected, delay time be- 
ing effective in smoothening out the shock-wave nature of 
the traveling wave. We also establish numerically that for 
the reaction-diffusion system the strength of the reaction 
term must not exceed a critical limit to allow travelling 
wave front solutions to exist for appreciable memory or 
relaxation effect. In view of the fact that the studies 
on reaction-diffusion and convection-diffusion with finite 
memory transport have been applied to forest fire |2~]] ] 
and population growth models U$ , Neolithic transitions 
p2[ and in several other areas under various approximate 
schemes @ 0, || || £| H |§ gj, (§, we be- 
lieve that the present exact solutions for the generalized 
Fisher and Burgers' problem are very much pertinent in 
this context. 



2[j+Jj^ = (U- Ul )(U-u 2 ) (40) 



to integrate to obtain finally 
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